function [u]= exsol(x,y);
% date le coordinate globali X=(x,y) 
% ...calcola la soluzione esatta 

% -----------------
% PROBLEMA STAZIONARIO (alfa=1,gamma=0)
% dominio quadrato [0,1]^2
% u = sin(x*pi).*sin(y*pi); 
% -----------------

% -----------------
% PROBLEMA STAZIONARIO (alfa=1,gamma=0)
% su dominio circolare centrato nell'origine, raggio 1
% u = 1 - (x.^2 + y.^2); 
% -----------------

% -----------------
% PROBLEMA STAZIONARIO (alfa=1,gamma=0)
% su dominio (poligonale) a forma di boomerang, 
% vedi dentro cartella "Meshes/saved_meshes"
% Caricare il file boomerang_mesh.mat
%
u = (y-x).*(y-x/2).*(y+x/2-1).*(y+x-2); 

% seguono gradiente-x, gradiente-y, laplaciano:

% -((x - 1)*(- 2*x^2 + 4*x + 5*y^2 - 6*y))/2
% -(5*x^2*y)/2 + (3*x^2)/2 + 5*x*y - 3*x + 4*y^3 - 9*y^2 + 4*y
% x^2/2 - x + (19*y^2)/2 - 15*y + 6
% -----------------
